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Dynamical chiral symmetry breaking (DCSB) is one of the keystones of low-energy hadronic phenomena. 
Dyson-Schwinger equations provide a model-independent quark-level understanding and correlate that with the 
behaviour of the pion's Bethe-Salpeter amplitude. This amplitude is a core element in the calculation of pion 
observables and combined with the dressed-quark Schwinger function required by DCSB it yields a valence-quark 
distribution function for the pion that behaves as (1 — x)"^ for s ~ 1, in accordance with perturbative analyses. 
This behaviour can be verified at contemporary experimental facilities. 



1. DCSB and the Gap Equation 

Dynamical chiral symmetry breaking (DCSB) 
is a signature characteristic of the strong inter- 
action spectrum and it is a feature of QCD that 
can only be understood via nonperturbative anal- 
ysis."'^ In this connection one particularly insight- 
ful tool is the QCD gap equation; i.e., the Dyson- 
Schwinger equation (DSE) [Q for the dressed- 
quark self-energy: 



S(p) ^ = 2^2 (17 • p -I- mbaro) 

+ Zi / g^D^,{p^q)^-f^S{q)Tl{q,p).{l) 



In Eq. (|l|): D^,y{k) is the renormahsed dressed- 
gluon propagator; T'^{q;p) is the renormahsed 
dressed-quark-gluon vertex; mbarc is the A- 
dependent current-quark bare mass that appears 
in the Lagrangian; and :— d^q/ rep- 
resents mnemonically a translationally-invariant 
regularisation of the integral, with A the reg- 
ularisation mass-scale. Also, Zi(C^,A^) and 
.^2(C^j^^) 8.re the quark-gluon- vertex and quark 
wave function renormalisation constants, which 
depend on the renormalisation point, C, and the 
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regularisation mass-scale, as does the mass renor- 
malisation constant: 



z,„(c^A2)=.z4(c^A2)/z2(c^A2), 



(2) 



with the renormahsed mass given by ra[C) 

^barc 

(A)/Z„,(C2,A2). 
The solution of Eq. (pi) has the form 



i^-pA{p\e)+B{p\C^) (3) 
^^-^[^7•p + M(p^C^)] , (4) 



where the functions A(p^,C^), B{p'^X^) express 
the effects of dressing induced by the quark's in- 
teraction with its own gluon field. Equation (|^) 
must be solved subject to a renormalisation con- 
dition, and in QCD it is practical to impose the 
requirement that at a large spacelike 



(5) 



A weak coupling expansion of the DSEs repro- 
duces every diagram in perturbation theory and 
in connection with Eq. (|l]) this means that at 
large spacelike-p^ the solution for massive quarks 
is, in Landau gauge. 



M(p2) 



I In 



pVA2 



QCD 



(6) 



^ These observations are also true of confinement but 
model-independent statements harder to make. 



2 



7m = 12/(33 — 2Nf), where rh is the 
renormaUsation-point-independent current-quark 
mass, and the mass renormahsation constant. 
Eq. (I), is Z,„(C2,A2) = [a(A2)/a(C2)]^-. At 
one-loop Zi = \ = Z{p^). In perturbation theory 
each contribution to the dressed-quark mass func- 
tion is proportional to the current-quark mass 
and hence M(p^) = in the chiral hmit, which is 
unambiguously defined by m = ||^. 

Dynamical chiral symmetry breaking is the ap- 
pearance of a C^) ^ solution of Eq. in 
the chiral limit. This guarantees a nonzero value 
of the vacuum quark condensate : 

~ (g-<z)° = Z4(C^A2)7V^trD T 5^9,0, (7) 

Jq 

where tTD identifies a trace over Dirac indices 
only and the superscript "0" indicates the quan- 
tity was calculated in the chiral limit. (NB. The 
factor of Zi here guarantees that the condensate 
is gauge-parameter- and cutoff-independent. Its 
omission yields a formula that is incorrect.) 

As I have just observed, B{p^,Q^) ^ is im- 
possible in perturbation theory. Hence a non- 
perturbative analysis of QCD's gap equation is 
required to explore this possibility. To arrive 
at model-independent conclusions a systematic, 
symmetry-preserving truncation scheme for the 
71-point functions in the gap equation must be 
used. 

1.1. Goldstone's Theorem 

At least one such scheme is known and it 
has been used to good effect. For example, the 
pion's Bethe-Salpeter amplitude has the general 
form 

ri(fc; P) = T^75 [iE^{k- P) 7 ■ m(fc; P)+ 

+ j-kk- PG^{k;P)+a^,k^P,H^{k;P)], 

(8) 

where P is the pion's total momentum and k is 
the relative momentum between the bound state's 
constituents. The scheme of Ref. ^ guaran- 
tees that the axial-vector Ward-Takahashi iden- 
tity (here written for m = 0, fc± = fc ± P/2): 

-iP^Ti^{k-P) 

= 5-i(fc+)75y +75y^-'(fc-), (9) 



where 

7 • kk^Gnik- P) - a^,k,HR{k- P)] + [. . .], 

(10) 

with FR{k;P), GR{k;P), HR{k-P) regular at 
P^ = 0, is satisfied at every order. This identity, 
which is fundamental to the successful application 
of chiral perturbation theory, leads to the follow- 
ing quark- level Goldberger-Treiman relations O] , 



which are exact in QCD: 








/X(fc;0) = 




(11) 


FR{k;Q)~ 


f2/OF,(fc;0) = 


A{e), 


(12) 


Gi?(fc;0)H 


-2/°G,(fc;0) = 


2A'(fc2) , 


(13) 


-ffi?(fc;0)H 


-2/OG,(fc;0) - 


0, 


(14) 



where is the pion's chiral-limit leptonic decay 
constant, obtained in general from 



Zi 



tr 



y757M5('z+)ri('z;^)5('z-) 



(15) 



with the factor of Z2 on the right-hand-side guar- 
anteeing that /tt is independent of the renormal- 
isation point, cutoff and gauge parameter. (A 
formula without Zi is incorrect.) Equations (pj])- 
( p^ ) form an essential part of the proof |^ that 
if, and only if, a theory exhibits DCSB then 
the homogeneous, isovector, pseudoscalar Bethe- 
Salpeter equation has a massless, P^ = 0, solu- 
tion and the axial-vector vertex is dominated by 
the pion pole for ~ 0; i.e., a proof of Gold- 
stone's theorem in QCD. 

1.2. A Pseudoscalar Meson Mass Formula 

Similar analysis of the m 7^ axial-vector 
Ward-Takahashi identity yields the following 
mass formula for flavour nonsinglet pseudoscalar 
mesons: 



(16) 



where: M.^ = 4 
current-quark masses of the meson's constituents; 
fn is the meson's leptonic decay constant, ob- 
tained via obvious analogy with Eq. ( p^ ; and the 



is the sum of the 
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in-meson condensate is defined via 



JH 

Z4 / tr 

Jq 



(17) 



J 



with, e.g., T"^ = (A^ + «A2)/2, where {A^' 
1, . . . , 8} are the Gell-Mann matrices, and (•)' de- 
noting matrix transpose. 

A number of observations are now in order. 



Equation (11) makes piain that in the presence 
of DCSB the magnitude of is set by the 
constituent-quark mass: M(p^ ~ 0). The same is 
true of {qq)f , which, in the chiral limit, is iden- 
tical to the vacuum quark condensate, Eq. (^. 
This hints at a corollary of Eq. ([l6|): in the chi- 
ral limit it yields the Gell-Mann-Oakes-Renner 
relation, a primary element in the application of 
chiral perturbation theory to light-quark mesons. 
Furthermore, while it is not readily apparent 
from the material presented heretofore, the DSE 
derivation of Eq. ([l6| ) assumes nothing about the 
magnitude of m^^ . Hence Eq. (|l^) is a single 
mass formula that unifies the light- and heavy- 
quark domains and it has consequently been used 
to prove that the mass of an heavy pseudoscalar 
meson rises linearly with the mass of its heaviest 
constituent H,^. 

2. Realising DCSB 

The truncation scheme of Ref . |^ also provides 
the foundation for a phenomenological model that 
has been used to very good effect in describing 
light-quark mesons . The model is founded 

on the assumption that the leading order term 
in the truncation of QCD's gap equation is well 
represented by a renormalisation group improved 
rainbow approximation p|: 



Sip) ^ = Z2 (17 -p + mbarc) 



A-^ 



A-^ 



+ J Q{{p - lY) - <D^i^.s{q)-i. , 

(18) 

where D^^''(fc) is the free gauge boson propaga- 
tor. In Eq. (ITsI), G{k'^) is an "effective coupling:" 



e(fc2) = 47ra(fc2) for P > 1 GeV^, where a{P) is 
the strong running coupling constant, but the be- 
haviour of C/(fc^) for fc^ < 1 GeV^; i.e., at infrared 
length-scales (> 0.2 fm), is currently unknown. 

It is immediately obvious from Eq. (|l|) that 
the effective coupling is a mnemonic for the con- 
tracted product of the dressed-gluon propagator 
and dressed-quark-gluon vertex. Hence the in- 
frared strength of the kernel in the gap equation 
can be characterised by an interaction-tension: 



1 

47r ./A 



'pQCD 



dkH' [A(fc2)_A(A2QCD)] (19) 



QCD 
2\ 



with A(fc2) = g{P)/k^, and ApQCD = IOAqcd 
the boundary above which perturbation theory 
is unquestionably valid. Studies show JTot that 
DGSB cannot occur for < 0.5 GeV^ - 9A^cd: 
which means that the existence of a dynamically 
generated, B ^ solution of Eq. (|l|) is impossi- 
ble without a significant infrared enhancement of 
the effective interaction. Reproducing observable 
phenomena requires 

>4GeV^~70A|cD; (20) 

i.e., a ten- fold enhancement over the critical 
value. The origin of this enhancement, whether 
in the gluon or ghost vacuum polarisation, or else- 
where, is currently unknown but is actively being 
sought (see, e.g., Refs. 0-|ll). That is a key to 
understanding the source of confinement as op- 
posed to simply representing it. 

References [pl|7|-p| employ a one-parameter 
model for g{k^, with that parameter, and the 
u = d- and s-current-quark masses, tuned to fit 
the experimental values of m^, /tt, mriK- This fit 
yields: the value of a'^ 



70A|cD in Eq. I 



iiGoV ^ 5.5 MeV, 



to]°°^ = 124 MeV; (21) 

and a value of the vacuum quark condensate: 

-(99>?GcV = (0-242 GeV)3, (22) 

for which a comparison with the calculated value 

-(99>TGeV = (0-245 GeV)3 (23) 

emphasises just how close one is to the chiral limit 
when dealing with pion observables. In addi- 
tion, the essential character of the drcssed-quark 
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propagator predicted by the model, which under- 
hes these condensate values, has recently been 
confirmed in numerical simulations of lattice- 
QCD 

Fixing the single parameter at this fitted value, 
many observables are predicted and the model 
achieves a r.m.s. error over predicted quantities 
of < 4% 1 15 . Furthermore, it is the only model 
to predict ||] a behaviour for the pion's electro- 
magnetic form factor that agrees with the results 
of a recent Hall C experiment The large- 

behaviour of the form factor can be obtained 
algebraically and one finds Q Q^F-^{Q'^) = 
const., up to logarithmic corrections, in agree- 
ment with the perturbative-QCD expectation. 
This result relies on the presence of pseudovec- 
tor components in the pion's Bethe-Salpeter am- 
plitude, which is guaranteed by the quark-level 
Goldberger-Treiman relations in Eqs. (|l2|), (p^. 

The unification of light- and heavy-meson 
masses via the mass formula in Eq. (jlj) has also 
been quantitatively explored using the model of 
Refs. [|[||-§. This is illustrated in Fig. || wherein 
the calculated mass of a uq pseudoscalar meson 
is plotted as a function of m^, with m^, fixed to 
a value corresponding to that in Eq. The 
fitted curve is, in MeV 



6000 



TUH = 83 + 500\/A^ + 310A', X = m^/AQCD,(24) 

with the renormalisation mass-scale C = 19 GeV 
and Aqcd = 0.234 GeV. 

In this figure the curvature appears slight but 
that is misleading: the nonlinear term in Eq. ( p^ ) 
accounts for almost all of m^r (the Gell-Mann- 
Oakes-Renner relation is nearly exact for the 
pion) and 80% of uik- NB. The dashed line in 
Fig. |l| fits the K, D, B subset of the data ex- 
actly. It is drawn to illustrate how easily one can 
be misled: without careful calculation one might 
infer from this apparent agreement that the large- 
rUq limit of Eq. (^6|) is already manifest at the s- 
quark mass. However, in reality, the linear term 
only becomes dominant for niq > 1 GeV, provid- 
ing 50 % of and 67 % of to^. The model pre- 



dicts 



,lGeV 



1.1 GeV and 



,lGoV _ 



4.2 GeV, 



values that are typical of Poincare covariant treat- 
ments of heavy-meson systems Q . 

A similar analysis of pseudoscalar mesons with 



4000 



> 



2000 




(MeV) 



Figure 1. Solid line: pseudoscalar uq me- 
son's mass as a function of m^, ^ — 19 GeV, 
with a fixed value of that corresponds to 
TO^^°^ = 5.5 MeV, Eq. The experimen- 

tal data points are taken from Ref. |18| as are 
the errors assigned to the associated heavy-quark 
masses. The dashed curve is a straight line drawn 
through the K, D, B masses. 

equally massive constituents has also been per- 
formed p3[ and this predicts |19{: 



mH„ 



2.2. 



(25) 



in agreement with a result obtained in recent 
quenched lattice simulations |2^] . The model cal- 
culation provides an understanding of the lattice 
result. It shows that the persistent dominance 
by the term nonlinear in the current-quark mass 
owes itself to a large value of the in-meson conden- 
sates for light-quark mesons; e.g., — (99)igoV ~ 
(0.32 GeV)'^ and thereby provides a confirma- 
tion of the large condensate version of chiral per- 
turbation theory. 



5 




Figure 2. DSE result for ||21|1 evolved 

to — 16 GeV^ using the first-order, nonsinglet 
renormalisation group equation, for direct com- 
parison with ttN Drell-Yan data |22) , filled circles. 
The filled squares illustrate the anticipated errors 
achievable in a proposed JLAB experiment 1 23 1 . 



3. Pion's Valence Quark Distribution 

I have illustrated that the DSEs furnish a sound 
theoretical description of pion properties, provid- 
ing a model-independent explanation of its essen- 
tial dichotomous nature as both a Goldstone bo- 
son and a low-mass bound state of massive con- 
stituents, and also that they supply an efhcacious 
phenomenological tool. This makes them ideal 
for exploring the nature of the parton distribu- 
tion functions in the pion. These functions, while 
they provide a direct measure of the pion's quark- 
gluon substructure via perturbative QCD, cannot 
be calculated using perturbative methods. 

The valence-quark distribution function. 
My (a;), is of particular interest because its point- 
wise behaviour is affected by aspects of con- 
finement dynamics, e.g., by the mechanisms re- 
sponsible for the finite extent and essentially 
nonpointlike nature of the pion. (NB. = d'^ 
in the fj-parity symmetric limit of QCD.) The 
low moments of these functions can be inferred 
from numerical simulations of lattice-QCD but 
such moments alone are insufficient to determine 



the functions: they are only weakly sensitive to 
the behaviour of the distribution functions on the 
valence-quark domain, x > 0.5. 

The DSE model used in illustrating the frame- 
work's ability to unify the small- and large-Q^ 
behaviour of the pion's form factors was used to 
calculate Uy (x) . This calculation of the ap- 
propriate "handbag diagrams" indicates that at 
a resolving scale go — 0.54 GeV= 1/(0.37 fm) va- 
lence quarks, with a mass of 0.3 GeV, carry 71 % 
of the pion's momentum, and yields the distribu- 
tion function depicted in Fig. ^. 

Over the entire range of go considered, the cal- 
culated distribution function is precisely fitted by 
a MSR form 

u^{x) = A„ a;"!-! (1 - o;)"^ (1 - e„yi + 7„a;),(26) 
with exemplary, calculated parameter values: 



qo (GeV) 




Vi 


V2 


e« 


7« 


0.54 


11.24 


1.43 


1.90 


2.44 


2.54 


2.0 


4.25 


0.97 


2.43 


1.82 


2.46 


4.05 


3.56 


0.89 


2.61 


1.62 


2.30 



(27) 



At qo ~ 2.0 GeV two low moments of the distri- 
bution function arc pl| : 









Gale. 


0.098 


0.049 


Exp. 


0.10 ±0.01 


0.058 ±0.004 


Latt. 


0.11 ±0.3 


0.048 ± 0.020 



(28) 



with the experimental results from Ref. |^ and 
the lattice results from Ref. (24j. NB. Given the 
evident disagreement between the DSE calcula- 
tion and the data, the agreement between the cal- 
culated moments emphasises the insensitivity of 
low moments to the x-dependence of the distri- 
bution function on the valence-quark domain. 

A material feature of the DSE result is the 
value of 772 — 2 because while perturbative QCD 
cannot be used to obtain the pointwise depen- 
dence of the distribution functions it does give a 
prediction for the power-law dependence at a; ~ 1. 
That prediction is 123]: 



pQCD: <(x) a {I - 



(29) 



as 



in agreement with the DSE result. However 
will have been anticipated from Fig. H this pre 
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n 
► 



Figure 3. Deep inelastic scattering from the pro- 
ton's TT-cloud could provide a means of measuring 
u^rix) at JLAB ||. 

diction disagrees with the extant experimental 
data p2[, an analysis of which yields: 



DrellYan: (x (1 - x) . 



(30) 



The disagreement is very disturbing because a 
verification of this experimental result would 
present a profound threat to QCD, even challeng- 
ing the assumed vector-exchange nature of the 
force underlying the strong interaction. 

The DSE study has refocused attention on 
this disagreement, and is the catalyst for a resur- 
gence of interest in Uy{x) and proposals for its 
remeasurement. One proposal that could use ex- 
isting facilities would employ the (Sullivan) pro- 
cess depicted in Fig. || at JLAB |23j, with the 
anticipated accuracy illustrated in Fig. |^. This 
process could also be used efficaciously at a fu- 
ture electron-proton collider to accurately probe 
Uy{x) on the valence-quark domain, as empha- 
sised by Fig. ||[|6|. 

4. Epilogue 

Dynamical chiral symmetry breaking (DCSB) 
is a keystone of hadron physics. It is the ef- 
fect responsible for turning perturbative current- 
quark masses into nonperturbative constituent- 
quark masses, and for ensuring that the pion 
is light while the pion's electroweak decay con- 
stant nevertheless sets a large mass-scale: 47r/7r ~ 
1 GeV. It is thereby the foundation for the suc- 
cessful application of chiral perturbation theory 
to low-energy hadronic phenomena. 

The QCD gap equation supplies a quark-level 
explanation of DCSB and, as one of the tower of 
Dyson-Schwinger equations (DSEs), unifies that 




0.85 
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Figure 4. Simulated errors for DIS events 
obtained in collisions of a 5 GeV electron beam 
and a 25 GeV proton beam with a luminosity 
of 10'^^ cm^^ s^^ and 10^ s of running. (Figure 
adapted from Ref. pGf.) 

with a Poincare covariant understanding of the 
structure of QCD's bound states and their inter- 
actions. In this way one finds, for example, that 
DCSB is also the main reason for the large tt — p 
and p — ai mass-splittings. 

The fact that perturbation theory is recovered 
via a weak coupling expansion of the DSEs pro- 
vides a tight constraint on the ultraviolet be- 
haviour of the calculated Schwinger functions 
that describe hadronic interactions. Therefore 
their use in predicting and correlating observ- 
ables, which necessarily sees the introduction of 
some model dependence, provides a means by 
which experimental results can be used to probe 
the infrared (long-range) behaviour of the quark- 
quark interaction; i.e., of exploring the mecha- 
nism of confinement. In this way the framework's 
adaptability to modelling is a material asset. 

As reviewed in Refs. [^^, modern applica- 
tions of the DSEs have met with substantial suc- 
cess, even at nonzero temperature and chemical 
potential, but challenges remain, of course. The 
gap equation makes clear that the existence of 
DCSB signals a significant enhancement of the 
strong coupling over the perturbative expectation 
on the infrared domain: fc^ < IGeV^. However, 
the mechanism in QCD that supplies that en- 
hancement is yet to be conclusively identified, and 
is being explored using DSE and lattice methods. 

A direct bound state treatment of the scalar 
meson sector is also wanting. However, at least 
one now understands why the lowest order non- 
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perturbative truncation of the kernels in the rele- 
vant integral equations (rainbow-ladder), so suc- 
cessful for pseudoscalar and vector mesons, fails 
for the scalars |2^. Improvements, systematic 
and/or imaginative, are being explored. 

This, after all, is a core issue in DSE studies 
and the primary point of criticism: the system 
of integral equations must be truncated and how 
does one judge, a priori, the fidelity of a given 
procedure? Addressing this open question is a key 
focus of contemporary research but the question's 
existence does not itself diminish the efficacy nor 
the value of modern DSE phenomenology. 
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